Accepted Version 1 



Thermal Model Calibration for Minor Planets Observed with 

WISE/NEOWISE 

A. Mainzer\ T. Grav^, J. Masiero^ 

J. Baueri'3, E. Wright^ R. M. Cutri^, R. S. McMillan^ M. Cohen^M. Ressler^R 

Eisenhardt^ 



ABSTRACT 



With the Wide-field Infrared Survey Explorer (WISE; Wright et aL 2010), 
we have observed over 157,000 minor planets (Mainzer et al. 2011). Included 
in these are a number of near-Earth objects. Main Belt Asteroids, and irregular 
satellites which have well-measured physical properties (via radar studies and in 
situ imaging) such as diameters. We have used these objects to validate mod- 
els of thermal emission and reflected sunlight using the WISE measurements, as 
well as the color corrections derived in Wright et al. (2010) for the four WISE 
bandpasses as a function of effective temperature. We have used 50 objects with 
diameters measured by radar or in situ imaging to characterize the systematic 
errors implicit in using the WISE data with a faceted spherical NEATM model 
to compute diameters and albedos. By using the previously measured diameters 
and H magnitudes with a spherical NEATM model, we compute the predicted 
fluxes (after applying the color corrections given in [Wright et aL]|2010D in each of 
the four WISE bands and compare them to the measured magnitudes. We find 
minimum systematic flux errors of 5-10%, and hence minimum relative diameter 
and albedo errors of ~10% and ~20%, respectively. Additionally, visible albedos 
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for the objects are computed and compared to the albedos at 3.4 /im and 4.6 
/im, which contain a combination of reflected sunlight and thermal emission for 
most minor planets observed by WISE. Finally, we derive a linear relationship 
between subsolar temperature and effective temperature, which allows the color 
corrections given in Wright et al. (2010) to be used for minor planets by comput- 
ing only subsolar temperature instead of a faceted thermophysical model. The 
thermal models derived in this paper are not intended to supplant previous mea- 
surements made using radar or spacecraft imaging; rather, we have used them to 
characterize the errors that should be expected when computing diameters and 
albedos of minor planets observed by WISE using a spherical NEATM model. 



1. Introduction 



The Wide-field Infrared Survey Explorer (WISE) is a NASA Medium-class Explorer 
mission designed to survey the entire sky in four infrared wavelengths, 3.4, 4.6, 12 and 22 



/im (denoted Wl, W2, W3, and WA respectively) (Wright et al. 2010 Liu et al. 2008 



Mainzer et al. 2005). The final mission data products are a multi-epoch image atlas and 
source catalogs that will serve as an important legacy for future research. While WISE's 
primary science goals are to find the most luminous galaxies in the entire universe and 
to find the closest and coolest stars, the survey has yielded observations of over 157,000 
minor planets, including Near-Earth Objects (NEOs), Main Belt Asteroids (MB As), comets, 
Hildas, Trojans, Centaurs, and scattered disk objects ( Mainzer et al.||2011[ ). This represents 
an improvement of nearly two orders of magnitude more objects observed than WISE's 
predecessor mission, the Infrared Astronomical Satellite (IRAS; Tedesco et al.|[l988 Matson, 



D. 1986). The WISE survey began on 14 January, 2010, and the mission exhausted its 
primary tank cryogen on 5 August, 2010. Exhaustion of the secondary tank and the start of 
the NEOWISE Post- Cryogenic Mission occurred on 1 October, 2010, and the survey ended 
on 31 January, 2011. 

Infrared observations of all classes of minor planets are useful for determining size and 
albedo distributions, as well as thermophysical properties such as thermal inertia, the mag- 



nitude of non-gravitational forces, and surface roughness (e.g., Tedesco et al. 2002 Trilling 



et al. 


2010 


Harris et al. 


2009) 



some have physical properties such as diameter, albedo, and subsolar temperature measured 
by independent means such as radar observations, in situ spacecraft visits, or stellar occul- 
tations. We can obtain diameters and albedos for the full sample of minor planets observed 
by WISE, but this requires characterization of the systematic errors associated with using 
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thermal models (c.f. Harris et al.||2011 ). In this work, we have used a set of objects with well- 
known, independently measured diameters and H magnitudes to test the ability of thermal 
models created using WISE data to accurately compute diameter and albedo. 

Figure [T] shows model spectral energy distributions for minor planets with different ef- 
fective temperatures with the WISE bandpasses overplotted. The color corrections given 



by Wright et al. (2010) allow the zero points to be corrected to account for the effect of 



flux generated by objects with varying effective temperature being observed through non- 
monochromatic bandpasses. The absolute zero points and isophotal wavelengths of the WISE 



passbands (Wright et al. 2010 Tokunaga & Vacca 2005) are calibrated using the particular 



Kurucz photospheric spectrum for Vega detailed by Cohen et al. ( 1992 ) and validated abso- 



lutely to 1.1% by Price et al. (2004). The width of the WISE passbands (particularly W?>) 



means that it is necessary to apply a color correction for sources with a different spectrum 
than Vega. For the vast majority of WISE extra-solar sources, these color correction terms 
are small (only a few percent for a K2V star in W?), for example). However, for objects such 
as minor planets with effective temperatures as low as 100 K, the flux correction factor for 
W2) exceeds a magnitude ( Wright et al.|2010 ), and readers are strongly encouraged to consult 
the WISE Explanatory Supplement and Wright et al. (2010) for precise values. Wright et 



al. (2010) also find a discrepancy between red and blue cahbrators in bands W?> and WA that 



require a -8% and -|-4% adjustment to the zero point magnitudes in each band, respectively. 



The color corrections derived in Wright et al. (2010) were produced by integrating over the 



system's relative spectral response throughput, which was measured during ground-based 



tests prior to launch (Latvakoski et al. 2010). 



We have compared the measured WISE magnitudes to magnitudes derived from spher- 
ical thermal models created for 117 objects, all of which have diameters and/or albedos 
measured via in situ spacecraft measurements, radar studies, or stellar occultations. By 
comparing the predicted magnitudes to the as-measured magnitudes reported by the WISE 



pipeline, we have verified that the color corrections derived in Wright et al. (2010) as a func- 



tion of an object's effective temperature correctly reproduce observed physical parameters. 
We have used these objects to set limits on the systematic errors for diameters and albedos 
derived from WISE observations of minor planets when using a spherical thermal model and 



the Wright et al. (2010) color corrections. Finally, we have determined an empirical rela- 



tionship between an object's subsolar temperature and its effective temperature, which is 



necessary for selecting the appropriate color correction for Solar System objects from Wright 
eraL|(pM0l). 



^http://wise2.ipac.caltech.edu/docs/release/prelini/expsup/ 
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2. Observations 

We have assembled a list of objects with well-measured diameters and H magnitudes 
that WISE observed during the fully cryogenic portion of its mission (Table 1). These 
objects were chosen because their physical measurements were obtained using methods 
largely independent of infrared model parameters, such as radar imaging, direct measure- 
ments from spacecraft visits or fiybys, or stellar occultations. Diameters and albedos derived 
from other infrared observatories such as IRAS, the Spitzer Space Telescope, and ground- 



based observations were not used to verify the color corrections of Wright et al. (2010) 
due to the desire to produce an independent calibration without reference to other infrared 
observers' model parameter assumptions. However, for some objects, such as those ob- 



served by Shepard et al. (2010), the radar diameter measurements were compared with 
IRAS diameters and were found to be consistent. Of the ~400 radar-observed asteroids 
{http://echo.jpl.nasa.gov/index.html), a dozen or so spacecraft targets, and hundreds of oc- 
cultation targets, we identified WISE observations for 117 objects. In order to reduce the 
potential difficulties associated with applying a spherical thermal model to non-spherical 
objects, we eliminated from further consideration all objects with peak-to-peak magnitude 
differences larger than 0.3 magnitudes in W3; this reduced the sample to the 50 objects 
listed in Table 1. 

The observations of these objects were retrieved by querying the Minor Planet Center's 
(MPC) observation files to look for all instances of individual WISE detections of the desired 
objects that were reported using the WISE Moving Object Processing System (WMOPS; 



Mainzer et al. 2011 ). The resulting set of position/time pairs were used as the basis of a query 
of WISE source detections in individual exposures (also known as "Level lb" images) using 
the Infrared Science Archive (IRSA). In order to ensure that only observations of the desired 
moving object were returned from the query, the search radius was restricted to 0.3 arcsec 
from the position listed in the MPC observation file. Additionally, since WISE collected a 
single exposure every 11 seconds and observes each part of the sky an average of 10 times, 
the modified Julian date was required to be within 2 seconds of the time specified by the 
MPC. Artifacts were minimized by setting the following fiags: ccJlags = or p and ph_qual 
= A, B, or C (this fiag indicates that the source is likely to have been a valid detection). 
Objects brighter than W3 = 4 and W4 = magnitudes were assumed to have fiux errors 
equivalent to 0.2 magnitudes due to changes to the shape of the point spread function as the 
objects became saturated, and a linear correction was applied to the W3 magnitudes in this 
brightness regime (the WISE Explanatory Supplement contains a more detailed explanation). 
Per the Explanatory Supplement, objects brighter than W3 = —2 and WA = —6 were not 
used. Each object had to be observed a minimum of three times in at least one WISE band, 
and it had to be detected at least 40% of the time when compared to the band with the 
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maximum number of detections (usually, though not always, W3). The WMOPS system is 
designed to reject inertially fixed objects such as stars and galaxies in bands W3 and Wi. 
Nonetheless, the individual images at all wavelengths were compared with WISE atlas coadd 
and daily coadd source lists to ensure that inertially-fixed sources such as stars and galaxies 
were not coincident with the moving object detections. This check is particularly important 
in bands Wl and W2 where the density of background objects (and hence the probability 
of a blended source) is higher than at longer wavelengths. Any remaining blended sources 
in bands Wl and W2 were removed. Some objects were observed at multiple epochs, and 
observations separated by more than three days were modeled separately. 



3. Thermal Model and Reflected Sunlight Fits 



The thermal fiux from an airless sphere is given by 

"17/2 r2iT 



sin 9 cos 9d9 



where A is the object-to-observer distance, e is the emissivity, R is the object's radius, B is 
the Planck function, 9 is the angle from the sub-observer point to a point on the asteroid such 
that 9 is equal to the solar phase angle a at the sub-solar point, and is an angle measured 
around the sub-observer point such that = at the sub-solar point (e.g, Bhattacharya et al. 
2010 Harris et al.|2009 Delbo fc Harris|2002 ). In order to compute this fiux, the temperature 
distribution across the body must be computed. In the Standard Thermal Model (STM) 



of Lebofsky & Spencer (1989), the temperature of an asteroid is assumed to be maximum 



at the subsolar point and zero on the point opposite of this; this is the case of an object 
with zero thermal inertia. In contrast, in the Fast Rotating Model (FRM) ( [Lebofsky et al. 



1978 Veeder et al. 1989; Lebofsky & Spencer 1989), the asteroid is assumed to be rotating 



much faster than its cooling time (i.e. high thermal inertia), resulting in a constant surface 
temperature across all longitudes. The so-called "beaming parameter" was introduced by 



Lebofsky (1986) in the STM to account for the enhancement of thermal radiation observed 



at small phase angles. The near-Earth asteroid thermal model (NEATM) of Harris (1998) 
also uses the beaming parameter rj to account for cases intermediate between the STM and 
FRM models, such that 



T{9, 



Tgs [max(0, cos 9 cos a + sin 9 sin a cos < 



11/4 



where a the solar phase angle, and the sub-solar temperature Tgg is given by: 

> H 0.25 



T 



go(l - A ) 
rjea 



(2) 



(3) 
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The emissivity, e, is assumed to be 0.9 for all wavelengths (c.f. Harris et al. 2009), A 



is the Bond albedo, 5*0 is the solar flux, a is the Stefan-Boltzmann constant, and rj is the 
beaming parameter. In the STM, t] is set to 0.756 to match the occultation diameters of 1 
Ceres and 2 Pallas, while in the FRM, t] is equal to tt. With NEATM, 77 is a free parameter 
that can be fit when two or more infrared bands are available (or with only one infrared 
band if diameter or albedo are known a priori as in this paper). 

Each object was modeled as a set of triangular facets covering a spherical surface with 



diameter equal to the ground-truth measurement (c.f. Kaasalainen et al. 2004). Model mag- 



nitudes were computed for each WISE measurement, ensuring that the correct Sun-observer- 
object distances were used for each one. The temperature for each facet was computed, and 



the Wright et al. (2010) color corrections were applied to each facet. The emitted thermal 



flux for each facet was calculated using NEATM along with the bandcenters and zero points 
given in Wright et al. (2010); the temperature at the anti-subsolar point was set to 3 K, so 



the facets closest to this point contribute little flux. 



The objects' absolute magnitudes (H) were taken from Warner, Harris & Pravec (2009) 



and Pravec et al. (2006) when available; otherwise, the values were taken from the MPCs 



orbital element files. Although Trilling et al. (2010), Parker et al. (2008) and Juric et al. 



(2002) have applied a 0.3 magnitude offset to the H values for objects in their respective 



works, we did not apply the offset to the H magnitudes ( Spahr|[2bll ). The offset was caused 
by the dominance of unfiltered CCD photometry in the MPC datafiles, largely from the 
LINEAR survey; as more filtered measurements have become available from other surveys 
in recent years, the offset is no longer appropriate to use. We will revisit these H values as 
improved measurements become available. 

Diameters and albedos computed from the combination of diameter and H from the 
radar, spacecraft, or occultation measurements are given in Table 1. In addition, the Satur- 



nian moon Phoebe had a measurement of subsolar temperature (Spencer et al. 2004) that 



was also used in the thermal model. Phoebe's H value was taken from Grav et al. (2003 ), and 



its G value from Bauer et al. (2006). If an object had no prior independent measurement of 



geometric albedo from an in situ measurement, it was computed using the relationship 



Pv 



1329 ■ 10-°-^" 
D 



(4) 



where H is the object's absolute magnitude and D is its diameter (Harris & Laggeros|[2002| . 
For objects with three valid measurements in two or more WISE bands dominated by thermal 
emission, 77 was determined using a least squares minimization. 



In general, minor planets detected by NEOWISE in bands Wl and W2 contain a mix of 
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reflected sunlight and thermal emission. Thus, it was necessary to incorporate an estimate 
of reflected sunlight into the thermal model in order to use data from bands Wl and W2. 
In order to compute the fraction of reflected sunlight in bands Wl and W2, it was also 
necessary to compute the ratio of the infrared albedo pm to the visible albedo pv- We 
make the simplifying assumption that the reflectivity is the same in both bands Wl and 
W2, such that pjr = ^3.4 = p^^', the vahdity of this assumption is discussed below. The 
geometric albedo pv is deflned as the ratio of the brightness of an object observed at zero 
phase angle (a) to that of a perfectly diffusing Lambertian disk of the same radius located 
at the same distance. The Bond albedo (A) is related to the visible geometric albedo pv by 
A ^ Ay = qpvi where q is the phase integral and is deflned such that q = 2 J ^{a)sin{a)da. 
$ is the phase curve, and q = 1 for $ = max {0, cos [a)). G is the slope parameter that 



describes the shape of the phase curve in the H — G model of Bowell et al. (1989) that 
describes the relationship between an asteroid's brightness and the solar phase angle. For 
G = 0.15, q = 0.384. Since q is never larger than about 2/3 for any observed value of G, 
the geometric albedo pv is not limited to <1.0; it is possible to have a value as large as 
~1.5 without violating conservation of energy. Conversely, values of pv approaching 1.0 still 
amount to integral reflectivity (Bond albedo) of ~60%, less than newly fallen snow, or other 
analogs of very high albedo. 

We make the assumption that pm obeys these same relationships, although it is possible 
that the phase integral described above varies with wavelength, so what we denote here as 
Pjfi for convenience may not be exactly analogous to pv- The flux from reflected sunlight was 



computed for each WISE band using the lAU phase curve correction (Bowell et al. 1989): 



=iJ + 51og(i?A)-2.51og[(l-G)<l'i(a) + G$2(a)], (5) 



where V is the predicted apparent magnitude; R is the heliocentric distance in AU; A 
is the geocentric distance in AU; a is the Sun-observer-object angle; G is slope parameter 
for the magnitude-phase relationship; and 



where 



= exp[-Aitan(-a)^'] 



i = l,2,Ai = 3.33, Bi = 0.63, A2 = 1.87, B2 = 1.22. 



(6) 
(7) 



Unless a direct measurement of G was available from Warner, Harris & Pravec (2009) 



or 



Parker et al. (2008), we assumed a G value of 0.15. Finally, the weighted averages of the 



model magnitude, measured WISE magnitude and Tgs were computed for each object by 
grouping together observations no more than a three day gap between them. This was done 
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to ensure that NEOs, which can have significant changes in distance over short times, were 
modeled accurately. 

Error bars on the model magnitudes and subsolar temperatures were determined for 
each object by running 100 Monte Carlo trials that varied the objects' H values by the error 
bars given in Table 1, their diameters by the error bars specified in Table 1, and the WISE 
magnitudes by their error bars using Gaussian probability distributions. The minimum 
magnitude error for all WISE measurements fainter than W?) = 4 and WA = 3 magnitudes 
was 0.03 magnitudes, per the in-band repeatability measured in Wright et al. (2010). For 
objects brighter than W3 = 4 and W4 = 0, the error bars were increased to 0.2 magnitudes, 
as these magnitudes represent the onset of saturation. Additionally, a linear correction was 
applied to objects with —2 < W3 < 4. The effect is thought to be caused by changes 
in the point spread function observed when images start to enter saturation (see the WISE 
Explanatory Supplement for a discussion of the effects of saturation on photometry). We have 
set the error bars to 0.2 magnitudes right at these limits in order to conservatively capture 
residual errors not fixed by the simple linear correction applied. Objects with W3 < —2 and 
W4 < —6 were too bright to be used. The error bar for each object's model magnitude was 
equal to the standard deviation of all the Monte Carlo trial values. If a measurement of H 



was not available from Warner, Harris & Pravec (2009) or Pravec et al. (2006), the error in 



H was taken to be 0.3 magnitudes (Spahr 2011). 



Figure [2] shows an example of an asteroid's WISE magnitudes as well as the resulting 



thermal model fit. Figures 3 abed show the difference between the average measured WISE 
magnitudes and the model magnitudes for bands Wl through WA using the color corrections 
given in Wright et al. (2010). In addition, we adjusted the W3 effective wavelength blueward 



by 4% from 11.5608/xm to 11.0984 /im, the WA effective wavelength redward by 2.5% from 
22.0883 /im to 22.6405 /xm, and we included the -8% and +4% offsets to the W3 and WA 
magnitude zeropoints (respectively) due to the red-blue calibrator discrepancy reported by 
Wright et al. (2010). The weighted means of the differences between observed and model 



magnitudes shown in Figure 3abcd are given in Table 2, and they are near zero for most 
objects. The apparent trend below zero for the objects with Tgs larger than ~300 K could 
be due to any of the following causes: 1) imperfect knowledge of the system relative spectral 
response curves, particularly in W4; 2) inaccuracies in either H values or diameters for 
some of the objects; 3) the temperature distribution may not precisely follow that given in 
Equation 2. This is less important for objects observed at low phase angles. The warmer 
objects observed by WISE tend to be NEOs observed at higher phase angles, which could 
lead to the subsolar point contributing significantly more fiux in the shorter wavelengths 
and the nightside of the object more fiux at longer wavelengths. Nevertheless, most of the 
predicted magnitudes are in good agreement with the observed magnitudes, indicating that 
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the procedure given in Wright et aL (2010) for color correction is adequate. 



Although many of the calibrator objects are known to be non-spherical (e.g. from radar 
shape modeling or optical lightcurves), the WISE observations generally consisted of ~10-12 
observations per object uniformly distributed over ~36 hours (Wright et al. 2010 Mainzer 



et al. 2005), so on average, a wide range of rotational phases was sampled. Although the 



variation in effective spherical diameter resulting from rotational effects tends to be averaged 
out, caution must be exercised when interpreting effective diameter results using spherical 
models for objects that are known to have large-amplitude lightcurve variations. Figure |4] 
shows the distribution of peak-to-peak W3 amplitudes for the entire sample of 117 objects. 
We have compared the mean differences between observed and model magnitudes for the 
entire sample of 117 objects to the sample of only those 50 objects with peak-to-peak W3 
amplitudes < 0.3 magnitudes, and there is no significant difference in the result given in 
Table 2. Nevertheless, to avoid any potential difficulties associated with applying spherical 
models to elongated objects, we excluded objects with peak-to-peak W3 amplitudes from 



further analysis. (Figure 3abcd shows the difference between observed and model magnitudes 
for only the 50 low-amplitude objects). For the objects in Table 1 for which pole orientations 



could be found in the lightcurve database of Warner, Harris & Pravec (2009), we computed 



the visible lightcurve amplitude at the aspect observed by WISE and found that for these 
objects, the apparent amplitudes were <~0.3 magnitudes. 



4. Systematic Diameter and Albedo Errors 

The offsets and errors given in Table 2 can be regarded as the minimum systematic 
errors in magnitude for minor planets observed by WISE/NEOWISE. Since diameter is pro- 
portional to the square root of the thermal flux (Equation 1), the minimum systematic diam- 
eter error due to uncertainties in the color correction is proportional to one-half the error in 
flux. These magnitude errors result in a minimum systematic error of ~5-10% for diameters 
derived from WISE data; they are of similar magnitude to the diameter uncertainties of most 
of the underlying radar and spacecraft measurements, which are ~10% (references are given 
in Table 1). Albedo is proportional to (Equation 4), and so it is linearly proportional to 
flux to first order. Therefore, minimum systematic errors on albedos computed from WISE 
observations are ~10-20%, subject to the assumption that spherical effective diameters can 
be computed for non-spherical shapes (future work will model the objects as non-spherical 
shapes and will use shape models and rotational information produced by lightcurve and 
radar studies). These should be regarded as minimum errors in cases of good signal-to- noise 
detections when the beaming parameter and the infrared albedo can be fit. It should also 
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be noted that these error estimates apply only to objects as distant as Saturn, as the most 
distant object we have considered herein is Saturn's moon Phoebe. Objects observed by 
WISE at greater distances (and therefore lower temperatures) may be subject to additional 
errors. 



5. Converting Subsolar Temperature to Effective Temperature 



The color corrections specified in Wright et al. ( |2010 ) are given as a function of input 



spectra that are assumed to blackbodies of different effective temperatures (Te//). In order 



to use the Wright et al. (2010) color corrections by computing only the subsolar temperature 



rather than a faceted thermophysical model for each object, we have computed the rela- 
tionship between Tss and (Teff). The total flux {F{Wn), where n = 1, ...,4) was computed 
using NEATM by applying a color correction to each facet on the sphere based on the facet's 
blackbody temperature. The equivalent flux at the center wavelength of each bandpass is 
computed, F{Xisophot)'i since this is a monochromatic flux, the color corrections are unity. 
We then derive an effective flux correction (/e//) which is given by 

ophot)- \^) 

The flnal step is to flnd the Tg// that gives the blackbody flux correction f{Teff) = feff- 

Figure [5] shows the relationship between Tgg and T^ff] a least-squares flt to these points 
yields the relationship 

Teff{K) = 0.95Tss{K)-3m. (9) 
We excluded objects with more than 20% reflected sunlight in a given band because their 



color corrections will be dominated by the small corrections needed for a G star (Wright et 



al. 2010). A single object can contribute up to four points on this plot, one for each band it 



is detected in with less than 20% reflected sunlight. Equation 9 provides a shortcut method 
to flux correct the WISE magnitudes in lieu of performing the facet-by-facet correction 
described above with negligible additional error. 



6. Albedo as a Function of Wavelength 



For objects that were detected according to the criteria described above in either Wl 
or W2, we computed the albedo at these wavelengths {pm, assuming psAfim = P4.6Atm) in 



addition to pv- Although Trilling et al. (2010) and Harris et al. (2009) assume that the 



albedo at 3.4 /im is 1.4 times higher than pv, this result is based on seven S-type Main Belt 
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asteroids observed at 3 fim by Rivkin et al. (1997). Harris et al. (2011) use observations 



of a number of NEOs observed by Spitzer and find Pm/pv consistent with 1.4, using the 



relationship between i] and a defined by Wolters et al. (2008). We have not used the Wolters 



et al. (2008) relationship to derive r] and instead have allowed it to vary. Since we know 



the diameter and can derive pv from Equation 4, we fit t] independently for each object in 
Table 1. Figure [6] shows the histogram of beaming parameters, and Figure [7] shows ?7 as a 
function of solar phase angle for the objects described in Table 1; however, we expect to 
significantly improve upon this result in a future work by using the general population of 
asteroids observed by WISE over a wide range of phase angles. 



Figure [8] shows Pm/pv for the objects detected in Wl and W2; as shown in Table 1, they 
represent a mix of NEOs, MBAs, and irregular satellites. The weighted average of the ratio 
of pm/pv is 1.27 ± 0.37. In computing pm/pv, we have made the simplifying assumption 
that G does not vary with wavelength, although it is known that G increases slightly, from 



0.15 to 0.20, when going from B to R bands (Bowell & Lumme 1979). However, in some 



cases, the variation in albedo with wavelength could be due to material and/or chemical 



properties (c.f. Phoebe; Clark et al. 2005). A future work comparing both the infrared and 
visible albedos with known taxonomic classifications is in progress. 



7. Conclusions 



The calibration methodology described in this work is not unique to the WISE infrared 
data; however, the uniquely large set of minor planet observations afforded by WISE/NEOWISE 
has permitted characterization of the systematic errors produced when applying spherical 
NEATM models to a number of objects with previously measured diameters that span a wide 
range of populations. We have studied NEOs, Main Belt asteroids, and irregular satellites. 
In particular, the selection of an effective temperature to use when applying a literature 
color correction can cause large changes in the resultant best-fit diameter unless the relation 
between Tgs and Te// is well- understood or a complete faceted model is employed. 



The color corrections described in Wright et al. (2010) have been used to produce an 



estimate of the minimum systematic errors in magnitudes for minor planets detected by 



WISE/NEOWISE (Mainzer et al. 2011). We have shown that the minimum diameter errors 
resulting from residual uncertainties in the color corrections and assumptions regarding G, 
H, and albedo are comparable to the uncertainties in the diameters measured by radar or in 
situ spacecraft imaging (i.e. ~10%); albedos can be determined to ~20% when good signal- 
to-noise multi-band WISE measurements and visible data are available. However, we note 
that objects more distant than the Trojan asteroids may be subject to increased systematic 
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errors, as the most distant object studied in this work was Saturn's moon Phoebe. Our 
model includes an estimate of reflected sunlight computed in all four bands, and we have 
computed albedos at 3.4 and 4.6 /xm as well as at visible wavelengths, which should yield 
interesting compositional insights when compared with spectroscopic taxonomic data. Areas 
for future improvement of thermal models include studying the effects of observing at high 
phase angles, modeling non-spherical shapes, and allowing both G and albedo to vary as 
a function of wavelength. All observations were processed with the Pass 1 (Preliminary) 
version of the WISE data processing pipeline; as the final version of the pipeline becomes 
available, which incorporates many improvements to instrumental calibration, we will revisit 
the thermal models for these objects. We have also derived a simple linear relationship 
between subsolar temperature and the effective temperature used in Wright et al. (2010), 
which should facilitate appropriate choices of color corrections for the WISE bandpasses. 
The thermal models computed herein are not intended to supplant the diameters measured 
by radar, in situ imaging or occultations, but rather they provide insight into the ability of 
spherical NEATM thermal models to accurately determine diameter and albedo when applied 
to the general population of WISE-observed minor planets for which these parameters have 
not been previously determined. This demonstrates that the WISE dataset offers a powerful 
new means of characterizing physical parameters of minor planets with great accuracy. 



8. Acknowledgments 



We thank L. Benner, M. Busch and M. Shepard for useful discussions and for provid- 
ing diameter information on a number of objects in advance of publication. This publi- 
cation makes use of data products from the Wide-field Infrared Survey Explorer, which is 
a joint project of the University of California, Los Angeles, and the Jet Propulsion Labo- 
ratory/California Institute of Technology, funded by the National Aeronautics and Space 
Administration. This publication also makes use of data products from NEOWISE, which is 
a project of the Jet Propulsion Laboratory/California Institute of Technology, funded by the 
Planetary Science Division of the National Aeronautics and Space Administration. We thank 
our referee, A. W. Harris of Pasadena, California, for constructive comments that materially 
improved this work. We gratefully acknowledge the extraordinary services specific to NEO- 
WISE contributed by the International Astronomical Union's Minor Planet Center, operated 
by the Harvard- Smithsonian Center for Astrophysics, and the Central Bureau for Astronom- 
ical Telegrams, operated by Harvard University. We also thank the worldwide community 
of dedicated amateur and professional astronomers devoted to minor planet follow-up obser- 
vations. This research has made use of the NASA/IPAC Infrared Science Archive, which is 



- 13 - 



operated by the Jet Propulsion Laboratory, California Institute of Technology, under con- 
tract with the National Aeronautics and Space Administration. This research has made use 
of NASA's Astrophysics Data System. 



- 14 - 



REFERENCES 

Bauer, J. M., Grav., T., Buratti, B., Hicks, M. 2006 Icarus 184, 181 
Benner, L., 2010 private communication 

Benner, L., Ostro, S., Magri, C, Nolan, M., Howell, E., Giorgini, J., Jurgens, R., Margot, 
J.-L., Taylor, P., Busch, M., Shepard, M. 2008 Icarus 198, 294 

Benner, L. A. M., M. C. Nolan, S. J. Ostro, J. D. Giorgini, C. Magri, and J. L. Margot. 2006 
Icarus 182, 474 

Benner, L., Ostro., S., Rosema, K., Giorgini, J., Choate, D., Jurgens, R., Rose, R., Slade, 
M., Thomas, M., Winkler, R., Yeomans, D. 1999 Icarus 137, 247 

Bhattacharya, B., Noriega- Crespo, A., Penprase, B., Meadows, V., Salvato, M., Aussel, H., 
Prayer, D., Ilbert, O., et al. 2010 ApJ, 720, 113 

Bowell, E., & Lumme, K. In Asteroids. 1979, University of Arizona Press, Tucson, Arizona, 
p 132-169. 

Bowell, E., in Asteroids II, 1989 

Brozovic, M., Benner, L., Magri, C., Ostro, S., Scheeres, D., Giorgini, J., Nolan, M., Margot, 
J.-L., Jurgens, R., Rose, R. 2010 Icarus, 208, 207 

Busch, M., Ostro, S., Benner, L., Giorgini, J., Jurgens, R., Rose, R., Magri, C., Pravec, P., 
Scheeres, D., Broschart, S. 2005 Icarus, 181, 145 

Busch, M., Giorgini, J., Ostro., S., Benner, L., Jurgens, R., Rose, R., Hicks, M., Pravec, P., 
Kusnirak, P., Ireland, M., Scheeres, D., Broschart, S., Magri, S., Nolan, M., Hine, A., 
Margo, J.-L. 2007 Icarus, 190, 608 

Busch, M., Ostro., S., Benner, L., Brozovic, M., Giorgini, J., Jao, J., Scheeres, D., Magri, C., 
Nolan, M., Howell, E., Taylor, P., Margot, J.-L., Brisken, W. 2011 astro-ph/1101.3794 

Clark, R., Brown, R., Jaumann, R., Cruikshank, D., Nelson, R., Buratti, B., McCord, 
T., Lunine, J., Baines, K., Bellucci, G., Bibring, J.-P., Capaccioni, P., Cerroni, P., 
Coradini, A., Formisano, V., Langevin, Y., Matson, D., Mennella, V., Nicholson, P., 
Sicardy, B., Sotin, C., Hoefen, T., Curchin, J., Hansen, G., Hibbits, K., Matz, K.-D. 
2005 Nature 435, 66 



Cohen, M., Walker, R. G., Barlow, M. J., Deacon, J. R. 1992, AJ, 104, 1650 



- 15 - 



Delbo, M., Harris, A., 2002, Meteoritics & Planetary Science 37, 1929 

Durech, J., Kaasalainen, M., Herald, D., Dunham, D., Timerson, B., Hanus, J., Frappa, E., 
Talbot, J., et al. 2011 Icarus in press 

Giorgini, J. D., , S. J. Ostro, L. A. M. Benner, P. W. Chodas, S. R. Chcslcy, R. S. Hudson, M. 

C. Nolan, A. R. Klemola, E. M. Standish, R. F. Jurgens, R. Rose, A. B. Chamberlin, 

D. K. Yeomans, J.-L. Margot 2002 Science 296, 132 

Grav, T., Holman, M., Gladman, B., Aksnes, K. 2003 Icarus 166, 33 

Harris, A., 1998 Icarus, 131, 291 

Harris, A., & Lagerros, J., in Asteroids HI, Bottke, W., Cellino, A., Paolicchi, P., Benzel, 
R., eds. 2002 University of Arizona Press, Tucson, p. 205-218 

Harris, A., Mueller, M., Lisse, C., Cheng, A. 2009 Icarus, 199, 86 

Harris, A., Mommert, M., Hora, J., Mueller, M., Trilling, D., Bhattacharya, B., Bottke W., 
Chesley, S., Delbo, M., Emery, J., Fazio, G., Mainzer, A, Penprase, B., Smith, H., 
Spahr, T., Stansberry J., Thomas, C. 2011 AJ 141, 75 

Juric, M., et al. 2002, AJ, 124, 1776 

Kaasalainen, M., Pravec, P., Krugly, Y., et al. 2004 Icarus, 167, 178 

Keller, H. U.; Barbieri, C; Koschny, D.; Lamy, P.; Rickman, H.; Rodrigo, R.; Sierks, H.; 
A'Hearn, M. F.; Angrilh, F.; Barucci, M. A.; et al. 2010 Science 5962, 190 

Latvakoski, H.; Cardon, J.; Larsen, M.; Elwell, J. 2010 Proc. SPIE, 7731, 25 

Lebofsky, L., Veeder, G., Lebofsky, M., Matson, D., 1978 Icarus, 35, 336 

Lebofsky, L. A.; Sykes, M. V.; Tedesco, E. F.; Veeder, G. J.; Matson, D. L.; Brown, R. H.; 
Gradie, J. C.; Feierberg, M. A.; Rudy, R. J. 1986 Icarus, 68, 239 

Lebofsky, L., & Spencer, J., Asteroids II, University of Arizona Press, 1989, p. 128-147 

Liu, F.; Cutri, R.; Greanias, G.; Duval, V.; Eiscnhardt, P.; Elwell, J.; Heinrichsen, 1.; 
Howard, J.; Irace, W.; Mainzer, A.; Razzaghi, A.; Royer, D.; Wright, E. L., 2008, 
SPIE, 7017, 16 

Mainzer, A.; Eiscnhardt, P.; Wright, E. L.; Liu, F.; Irace, W.; Heinrichsen, I.; Cutri, R.; 
Duval, V. 2005, SPIE, 5899, 262 



- 16 - 



Mainzer, A., J. Bauer, T. Grav, J. Masiero, R. McMillan, R. Cutri, J. Dailey, E. Wright, et 
al., 2011 ApJ731, 53 

Magri, C, Ostro, S., Rosema, K., Thomas, M., Mitchell, D., Campbell, D., Chandler, J., 
Shapiro, I., Giorgini, J., Yeomans, D. 1999 Icarus, 140, 379 

Magri, C, Nolan, M., Ostro, S., Giorgini, J., 2007, Icarus, 186, 126 

Mahapatra, R, Benner, L., Ostro, S., Jurgens, R., Giorgini, J., Yeomans, D., Chandler, J., 
Shapiro, I. 2002 Planetary and Space Science 50, 257 

Matson, D., ed. The IRAS Asteroid and Comet Survey, 1986, JPL D-3698 (Pasadena: JPL). 

Mitchell, D., Ostro, S., Rosema, K., Hudson, R., Campbell, D., Chandler, J., Shapiro, I., 
1996 Icarus, 118, 105 

Muinonen, K.; Piironen, J.; Shkuratov, Yu. G.; Ovcharenko, A.; Clark, B. E. Asteroids III, 
Bottke, W., Cellino, A., Paohcchi, P., Binzel, R., eds. University of Arizona Press, 
Tucson, p. 123-128 

Ostro, S., Campbell, D., Shapiro, I. 1983 AJ, 88, 5650 

Ostro, S., Campbell, D., Hine, A., Shapiro, I., Chandler, J., Werner, C, Rosema, K. 1990 
AJ, 99, 2012 

Ostro, S., Campbell, D., Chandler, J., Shapiro, I., Hine, A., Velez, R., Jurgens, R., Rosema, 
D., Winkler, R., Yeomans, D. 1991 AJ, 102, 1490 

Ostro, S., Nolan, M., Margot, J.-L., Magri, C, Harris, A., Giorgini, J. 2001 Icarus, 152, 201 

Parker, A., Ivezic, Z., Juric, M., Lupton, R., Sekora, M., Kowalski, A. 2008 Icarus, 198, 138 

Porco, C, et al. 2003 Science, 299, 1541 

Porco, P., et al. 2005 Science, 307, 1237 

Price, S. D, Paxson, C, Engelke, C, Murdock, T. L., 2004 AJ, 128, 889 

Pravec, P., Schcirich, P., Kusnirak, P., Sarounova, L., Mottola, S., Hahn, G., Brown, P., 
Esquerdo, G., et al. 2006 Icarus 181, 63 

Rivkin, A., Howell, E., Clark, B., Lebofsky, L., Britt, D. 1997 LPI 28, 1183 

Romanishin, W., & Tegler, S., 2005, Icarus, 179, 523 



-17- 



Shepard, M., Margot, J.-L., Magri, C, Nolan, M., Schlieder, J., Estes, B., Bus, S., Volquard- 
scn, E., Rivkin, A., Benner, L., Giorgini, J., Ostro, S., Busch, M. 2006 Icarus, 184, 
198 

Shepard, M., et al. 2010 Icarus, 208, 221 

Shevchenko, V., & Tedesco, E., 2006 Icarus, 184, 211 

Spahr, T., 2011 private communication 

Spencer, J., Pearl, J., Segura, M. 2004 AGUFM P43B, 03 

Taylor, P., Howell, E., Nolan, M., Benner, L., Brozovic, M., Giorgini, J., Vervack, R., 
Fernandez, Y., Magri, C., Mueller, M. 2010 DPS 42, 5704 

Tedesco, E., Matson, D., Veeder, G., Lebofsky, L., 1988, Comets to Cosmology: Proceedings 
of the Third IRAS Conference, Springer Berhn, Heidelberg, Vol. 297, p. 19-26 

Tedesco, E., Noah, P., Noah, M., Price, S. 2002 AJ, 123, 1056 

Tholen, D. J., Eds., Asteroid Absolute Magnitudes V12.0. EAR-A-5-DDR-ASTERMAG- 
V12.0. NASA Planetary Data System, 2009. 

Thomas, P., Veverka, J., Simonelli, D., Helfenstein, P., Carcich, B., Belton, M., Davies, M., 
Chapman, C, 1994, Icarus, 107, 23 

Thomas, P. C; Bclton, M. J. S.; Carcich, B.; Chapman, C. R.; Davies, M. E.; Sullivan, R.; 
Veverka, J. 1996, Icarus, 120, 20 

Timerson, B., Durech, J., Aguirre, S., Benner, L., Blacnhette, D., Breit, D., Campbell, S., 
Campbell, R., et al. 2009 MPBu 36, No. 3, p. 98-100 

Tokunaga, A. T. & Vacca,W. D. 2005, PASP, 117, 1459 

Trilling, D., et al. 2010, AJ, 140, 770 

Veeder, G., Banner, M. S.; Matson, D. L.; Tedesco, E. P.; Lebofsky, L. A.; Tokunaga, A. T. 
1989 AJ, 97, 1211 

Warner, B., Harris, A., Pravec, P. 2009 Icarus 202, 134 

WiUiams, G., personal communication, 2011 

Wolters, S., Green, S., McBride, N., Davies, J. 2008 Icarus 193, 535 



- 18 - 

Wright, E. L., et al. 2010 A J, 140, no. 6 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 19 - 



Table 1. Spherical NEATM models were created for 50 objects ranging from NEOs to 
irregular satellites in order to characterize the accuracy of diameter and albedo errors 
derived from NEOWISE data. The diameters and H values used to fit each object from the 
respective source data (either radar, spacecraft imaging, or occultation) are given. Objects 
that are hsted twice were observed at two epochs separated by more than three days. When 

observations were separated by more than three days, a separate calculation was made. 



Object 


D (km) 


H 


Vv 


PiR 


Tss{K) 


V 


Ref 


5 


115 ±12 


6.9 


0.25 ±0.05 


0.33 ±0.07 


232.0 ±3.0 


0.99 ±0.07 


a 


8 


140 ± 14 


6.4 


0.26 ±0.05 


0.44 ±0.06 


285.0 ±4.8 


0.79 ±0.03 


a 


13 


227 ± 38 


6.7 


0.07 ±0.01 


0.04 ±0.18 


252.2 ± 13.6 


0.89 ±0.21 


c 


22 


143 ± 14 


6.5 


0.17 ±0.06 


0.29 ± 0.06 


226.0 ±3.2 


1.11 ±0.09 


a 


27 


118 ± 17 


7.0 


0.20 ± 0.03 


0.46 ± 0.24 


236.9 ± 5.5 


1.08 ±0.10 


c 


31 


280 ± 23 


6.7 


0.05 ±0.01 


0.08 ±0.01 


217.3 ±2.6 


0.88 ±0.05 


c 


36 


103 ±1 


8.5 


0.07 ±0.01 


0.06 ±0.01 


283.5 ±0.8 


0.83 ±0.01 


c 


39 


163 ± 16 


6.1 


0.25 ±0.06 


0.48 ±0.10 


260.8 ± 4.2 


0.83 ±0.07 


a 


47 


138 ± 13 


7.8 


0.07 ±0.02 


0.06 ±0.01 


227.0 ±3.2 


1.13 ±0.07 


d 


53 


115 ±8 


8.8 


0.04 ±0.00 


0.04 ±0.01 


274.7 ±2.9 


1.06 ±0.05 


c 


54 


142 ± 14 


7.7 


0.05 ±0.01 


0.07 ±0.01 


280.9 ± 5.5 


0.84 ±0.07 


a 


83 


84 ±8 


8.7 


0.09 ±0.04 


0.13 ±0.03 


259.1 ±4.5 


0.92 ±0.07 


c 


85 


163 ± 16 


7.6 


0.06 ±0.02 


0.07 ±0.04 


224.8 ±3.8 


1.03 ±0.07 


a 


94 


187.5 ±27 


7.6 


0.05 ±0.02 


0.06 ±0.02 


219.9 ±4.9 


1.09 ±0.10 


d 


97 


83 ±6 


7.7 


0.21 ± 0.04 


0.29 ± 0.02 


273.0 ± 3.4 


0.74 ±0.05 


b 


105 


119 ± 11 


8.6 


0.05 ± 0.02 


0.03 ±0.01 


294.5 ±5.1 


0.90 ±0.07 


b 


114 


100 ±9 


8.3 


0.09 ±0.02 


0.14 ±0.04 


263.3 ± 4.5 


0.98 ±0.07 


c 


114 


100 ± 16 


8.3 


0.09 ±0.03 


0.13 ±0.05 


251.6 ± 7.0 


1.01 ±0.12 


c 


128 


188 ± 29 


7.5 


0.05 ±0.03 


0.05 ±0.02 


247.8 ±6.3 


0.96 ±0.10 


c 


135 


77 ±8 


8.2 


0.15 ±0.06 


0.25 ±0.05 


223.7 ±2.8 


1.20 ±0.08 


e 


139 


164 ± 19 


7.9 


0.04 ±0.01 


0.06 ±0.01 


244.2 ± 4.3 


0.93 ±0.07 


d 


145 


151 ±23 


8.1 


0.04 ±0.01 


0.04 ±0.01 


236.8 ± 5.3 


1.09 ±0.11 


d 


194 


169 ± 18 


7.7 


0.05 ±0.02 


0.06 ±0.02 


237.0 ±4.4 


0.90 ±0.07 


c 


198 


57 ±7 


8.3 


0.26 ±0.06 


0.37 ±0.07 


245.0 ±5.7 


0.89 ±0.11 


d 


208 


45 ±5 


9.3 


0.17 ±0.06 


0.29 ±0.06 


224.6 ±7.4 


1.16 ±0.18 


a 


208 


45 ±5 


9.3 


0.16 ±0.05 


0.29 ±0.06 


231.1 ±6.7 


1.06 ±0.14 


a 


211 


143 ± 13 


7.9 


0.06 ± 0.02 


0.07 ±0.01 


217.4 ±3.0 


0.94 ±0.06 


d 


230 


109 ± 16 


7.3 


0.17 ±0.02 


0.30 ±0.11 


251.7 ±6.5 


0.97 ±0.10 


c 


266 


109 ±7 


8.5 


0.06 ±0.01 


0.05 ±0.01 


252.2 ±2.1 


0.92 ± 0.04 


d 


308 


144 ± 13 


8.2 


0.05 ±0.01 


0.08 ±0.02 


238.8 ±3.1 


1.08 ±0.06 


e 


345 


99. ±9 


8.7 


0.06 ±0.01 


0.06 ±0.01 


276.7 ±4.4 


0.93 ±0.06 


e 


350 


99.5 ± 5 


8.4 


0.08 ±0.02 


0.08 ±0.01 


227.4 ± 1.5 


0.87 ±0.03 


e 


444 


163 ± 36 


7.8 


0.05 ± 0.04 


0.06 ± 0.06 


229.7 ± 10.2 


0.95 ±0.15 


d 
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Table 1 — Continued 



Ohiert 


D (km) 


H 




Pv 






T (K) 






Ref 


488 


150 ± 21 


7.8 





06 ±0 


01 


0.06 ±0.02 


219.5 ±4.4 





83 ±0.07 


d 


522 


84 ±9 


9.1 





06 ±0 


03 


0.09 ±0.02 


215.4 ±4.5 





84 ±0.07 


e 


566 


134 ± 15 


8.0 





06 ±0 


02 


0.09 ±0.02 


218.3 ±3.6 





81 ±0.06 


e 


654 


127 ± 13 


8.5 





04 ±0 


01 


0.05 ±0.01 


250.7 ±3.9 





95 ±0.07 


d 


704 


312 ± 17 


5.9 





08 ±0 


02 


0.09 ±0.01 


225.8 ±2.0 





88 ±0.03 


d 


704 


312 ±30 


5.9 





08 ±0 


03 


0.10 ±0.01 


224.4 ±2.9 





84 ±0.05 


d 


925 


58 ±6 


8.3 





25 ±0 


05 


0.38 ±0.10 


244.6 ±4.5 





90 ±0.08 


a 


951 


12 ±1 


11.5 





33 ±0 


04 


0.47 ±0.07 


268.4 ±4.6 





68 ±0.05 


g 


1512 


65 ±7 


9.6 





06 ±0 


03 


0.00 ±0.00 


243.4 ±4.7 





66 ± 0.06 


e 


1627 


9± 1 


12.9 





15 ±0 


03 


0.22 ±0.01 


256.8 ±4.4 


1 


23 ±0.10 


h 


1866 


8.7 ± 1 


12.7 





19 ±0 


07 


0.30 ±0.05 


209.4 ± 7.8 


1 


56 ±0.25 


i 


2867 


5.3 ± 1.2 


13.4 





28 ±0 


13 


0.56 ±0.29 


234.5 ± 23.8 


1 


33 ±0.55 


j 


7335 


0.9 ±0.4 


17.0 





31 ±0 


30 


0.40 ± 0.30 


276.4 ±43.0 


1 


40 ±0.98 


k 


68216 


1.4 ±0.2 


16.3 





27 ±0 


10 




300.6 ± 13.4 


1 


02 ±0.22 


1 


68216 


1.4 ±0.2 


16.3 





31 ±0 


19 




305.0 ±21.7 


2 


24 ±0.55 


1 


164121 


1.1 ±0.3 


16.2 





36 ±0 


20 


0.50 ±0.26 


344.1 ±43.2 





88 ±0.39 


m 


Himalia 


136.0 ± 12 


7.9 





07 ±0 


01 


0.07 ±0.01 


187.3 ±4.5 





88 ±0.10 


n 


2005 CR37 


1.0 ±0.2 


18.9 





05 ±0 


02 


0.06 ±0.03 


300.4 ±26.1 





98 ±0.31 


o 


2008 EV5 


0.4 ±0.0 


19.7 





13 ±0 


01 


0.18 ±0.01 


331.1 ±6.8 


1 


96 ±0.15 


P 


Phoebe 


213.2 ±1.3 


6.6 





09 ±0 


02 


0.06 ±0.00 


123.4 ±0.3 


1 


23 ±0.01 


q 
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Table 2. The mean offset (Am) and standard deviation (am) of the observed - model 
magnitudes for the cahbrator objects shown in Figure 3abcd[ N is the number of 

observations used. 



Band N Am crm 



Wl 46 0.00 0.02 

W2 47 0.01 0.10 

W3 52 -0.11 0.14 

W4 50 -0.03 0.17 




Wavelength (microns) 



Fig. 1. — : This figure illustrates the need for the color corrections given by Wright et al. 



(2010) in order to capture the shift in zero point and effective wavelength as a function 



of effective temperature: the WISE bandpasses are broad, particularly Wi. The WISE 
bandpasses are shown as shaded bars, and normalized fluxes are plotted for a range of 
different blackbody temperatures, illustrating the necessity for a color correction that varies 
as a function of an object's effective temperature. 
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Fig. 2. — : WISE observations as a function of modified Julian date as well as a spherical 
NEATM model fit are shown for a typical calibrator object, (208) Lacrimosa. NEOWISE 
detections of this asteroid span ~17 hours. The Wl points are shown in dark blue; W2 in 
cyan, W3 in green, and W4 in red. 
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Fig. 3abcd. — : With WISE, we have observed objects with diameters that have been mea- 
sured either with radar observations, in situ spacecraft visits, or by stellar occultations (see 



Table 1). After applying the color corrections specified in Wright et al. (2010), we can com- 



pare the observed WISE magnitudes to the fluxes predicted by a thermophysical model for 
bands Wl through Wi (Figures (a)-(d), respectively). The dashed line shows the weighted 
mean value of all the points; the dotted lines are the 1-a errors. Objects with WISE mea- 
surement errors that are less than 0.1 magnitudes are shown as red dots; objects with errors 
greater than 0.1 magnitudes are shown as black dots. 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

W3 Peak to Peak Amplitude Variation (mag) 

Fig. 4. — : The distribution of peak-to-peak amplitudes in W3 for 117 objects with in- 
dependently measured diameters peaks at ~0.25 magnitudes. Objects with peak-to-peak 
amplitudes larger than 0.3 magnitudes were excluded from our analysis. 
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Fig. 5. — : The correlation between subsolar temperature and effective temperature is well- 
described by a linear relationship. Wl — WA detections are shown in blue circles; W2 as 
cyan triangles; W3 as green squares; and W4 as red inverted triangles. Each object can 
contribute up to four points on the plot, depending on the number of bands in which it was 
detected. 
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Fig. 6. — : The beaming parameter rj values resulting from our NEATM fits are compared 
with solar phase angle a for the 50 objects listed in Table 1. The mean value for t] = 
0.96 ±0.28. 
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Fig. 

7] = —0.00054a + 0.97 and is plotted as a red dashed line. The relationship from Wolters et 
al. (2008) is given by ry = 0.013a + 0.91 and is shown as a blue dash-dot line. Future work 



7. — : The beaming parameter vs. phase angle. The best linear fit to our data is 



will examine the relationship between rj and a using the full WISE dataset, which includes 
many more objects over a wide range of phase angles. 



-27- 




1.0 1.5 2.0 
pIR / pV 



3.0 



Fig. 8. — : The ratio of the albedo at VFl and W2 (we assume the albedo is the same in 
both of these bands) compared to the albedo at visible wavelengths as a function of subsolar 
temperature. The dashed line indicates the weighted mean value oipm/pv = 1-27 ± 0.37. 



